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1 Introduction 



We study the existence of traveling wave solutions for lattice differential equations (LDEs) 
by means of a perturbation argument and Fredholm theory for mixed type functional 
differential equations, In particular, we prove persistence of traveling waves for a general 
class of lattice differential equations with bistable nonlinearity. Consider the following 
equation, 

f Uj = di(u j+ i - 2Uj + Uj-i) + d 2 (u j+2 - 2Uj + Uj- 2 ) ~ fa(Uj), j G Z 
\ fa( u ) = u ( u — a )( u — 1), O-i d>x, G? 2 G 1R- 

Our primary interest is in competing interactions between first and second nearest neigh- 
bors when d\ < and d 2 < 0. We develop a general technique for continuation of 
solutions of vector dissipative lattice differential equations and obtain results on existence 
of traveling front solutions for ( 11. II) when d\ < 0, < — d 2 <C 1 and when d 2 < and 
< 1. 

Our contribution is to develop techniques based upon the implicit function theorem that 
are applicable for vector equations that are similar to that developed by Bates, Chen, and 
Chmaj [1] for scalar equations. Whereas in [1] the limiting system is the traveling wave 
equation associated with the PDE u t = u xx — f(u), we consider, through the use of the 
Fredholm theory for mixed type functional differential equations [12] , limiting equations 
that may correspond to lattice differential equations. Among the chief motivations in this 
work (and in pQ) for the use of implicit function theorem based techniques is the desire 
to handle cases in which there does not exist a comparison principle. 

Traveling wave solutions to (II. ip have been extensively studied when d\ > and d 2 = 0. 
In particular, the work of Weinberger based upon the development of an abstract compar- 
ison principle is applicable to both PDEs and LDEs, although primarily for monostable as 
opposed to bistable problems. Zinner proved existence of traveling fronts using topological 
fixed point results [19] and stability [H] in the bistable case. A general stability theory 
was developed by Chow, Mallet- Paret, and Shen [5] and Shen employed comparison prin- 
ciple techniques to prove results on existence, uniqueness, and stability of traveling fronts 
in which / = f(u,t) may depend periodically on t. More recently Chen, Guo, and Wu 
developed a framework for existence, uniqueness, and stability of bistable equations in 
periodic media [4] and Hupkes and Sandstede [9] prove the existence of traveling pulse 
solutions for discrete in space Fitz-Hugh Nagumo equations that occur when coupling a 
relaxation variable to the discrete Nagumo equation ( 01.11) with d\ > and d 2 = 0). 

Associated with traveling waves for ( II. ip when d\ > and d 2 = is the mixed type 
functional differential equations 

-V(0 = difo>(£ - 1) - 2^(0 + <p{£ + 1)) - /(?(£)) 

which results from the traveling wave ansatz Uj(t) — (p(j — ct). Among the important 
contributions to the study of these types of equations is the pioneering work of Rustichini 
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[T5| [16], the development by Mallet-Paret of a Fredholm theory for linear mixed type 
function differential equations [T2] and its use to understand the global structure of trav- 
eling wave solutions [13]. Exponential dichotomies for these equations were investigated 
in [6] and [H] and center manifold theory and Lin's method were developed in [7] and [8], 
respectively. 

The case in which di < and d 2 = was investigated in [17] and [2]. In [T7j a model 
was developed for the dynamics of twinned microstructures that arise in martensitic phase 
transformation, e.g., in shape memory alloys, which led to (11 .ip in an overdamped limit. 
Subsequently, the bistable nonlinearity f{u) = u — H{u — a), H the Heaviside step func- 
tion, was employed and transform techniques were utilized to determine waveforms and 
wavespeeds. In [2] the cubic nonlinearity was employed and the problem was converted to 
a periodic media problem so that the results of pi] could be applied. A wealth of traveling 
wave solutions of both bistable and monostable type were revealed. Similar techniques 
may be used to determine traveling fronts when d 2 < and d\ = which results in two 
decoupled systems of equations. In [T7J [2] one of the essential ideas (see also [3]) was to 
convert to a system in terms of odd and even lattice sites. This effectively allows us to con- 
sider connecting orbit problems between vector equilibria as opposed to connecting orbit 
problems between time independent spatially periodic solutions. Existence and structure 
of traveling fronts for higher space dimension versions of (II. II) was recently investigated 
in [10] using comparison principle and continuation techniques. 

This paper is organized as follows. In section [2] we present some of the notation we will 
employ and background on Fredholm theory from [12] for linear mixed type functional 
differential equations. In addition, we summarize two approaches to the existence of 
traveling wave solutions in lattice differential equations. The first due to Chen, Guo, and 
Wu [1] provides existence, uniqueness and stability results for traveling front solutions 
of (II .ip when d\ and d 2 are positive. The second is due to Bates, Chen, and Chmaj 
[1] and provides existence of traveling front solutions when di + Ad 2 > 0. Section [3] 
contains our main results and establishes the persistence of traveling wave solutions for 
vector equations. In particular, we consider systems of lattice equations and allow, under 
certain non-restrictive conditions, general limiting systems. In section H] we consider the 
application of general results in section |3] to the existence of traveling fronts to (II .ip for 
values of d\, d 2 which even after rewriting as a system (equivalently in a periodic media) 
do not possess a comparison principle. 
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2 Preliminaries and Notation 



2.1 Fredholm Alternative for Lattice Differential Equations 

If X, Y are Banach spaces with norms || ■ \\x, \\-\\y respectively, then we let L(X, Y) denote 
the Banach space of bounded linear operators T : X — > Y . Denote the kernel and range 
of T G L(X, Y) by 

K(T) = {iG X\Tx = 0} and R(T) = {y G Y\y = Tx for some x G X}. 

Recall that T is a Fredholm operator if T satisfies the following: 

(i) K(T) is finite dimensional in X; 

(ii) R{T) is closed and also finite dimensional in Y. 
The Fredholm index of T is defined as 

ind(T) = dim(K(T)) — codim(R(T)) . 

II " ||l 2 ;|| ' an d || ■ Wh 1 denote the norms of the spaces L 2 (R, R N ), L°°(R, R N ) and 
H 1 (R,R N ), respectively Also, (f,g) := J[figi + ... + fuQ^dx, where / = (f u f N ), g = 
(gi,...,gN) e L 2 (R,R N ). We call / _L g if and only if (f,g) = 0. K L (T) denotes the 
orthogonal complement of the kernel K(T), that is, K L {T) = {/ G X : (/, g) — for all 
9 e K(T)}. 

In [12], Mallet-Paret investigated the Fredholm alternative for the following functional 
differential equations of mixed type, for 1 < p < oo, 

u\x) =^2A j (x)u(x + r j ) + f(x),ue W 1,P (R, I N ) (2.1) 
i=i 

where I is some bounded interval, r\ = 0, and rj ^ r^, 1 < j < k < Ni, Ni > 2. 

We may write it as . 

m'(x) = Lu + f(x), (2.2) 

and we have the homogeneous equation 

u\x) = Lu. (2.3) 

If Aj (x) is a constant matrix, which is independent of x, we denote it by Aj ;0 and then 
we may write equation (12. 3ft as 

u (x) = L u. (2.4) 
Define A Lo (s) = si - J^fli A jfl e sr i . We say (E3D is hyperbolic if A io (i0) ^0,6 eR. 
In the case, lim AAx) = Af for 1 < j < N t . Let L^u := ^J^i Af u(x + r,). We 

say equation (12.31) is asymptotically hyperbolic if equations (12. 4ft with replacing Lq by 
are hyperbolic. 

Define by . . , , 

A L u = -cu'(x) - Lu. (2.5) 

We recall Theorem A in Mallet- Par et's paper |12j: 
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Theorem 2.1. (See Jfflj) For each p with 1 < p < oo, Al is a Fredholm operator from 
W 1,p to LP provided that equation —cu'(x) = Lu is asymptotically hyperbolic. 

We note here that for linear mixed type functional differential equations the standard 
formula for computation of the Fredholm index is generally not valid, but this is remedied 
using the spectral flow formula (see [12] Theorem C). 

2.2 Traveling waves for Bistable Dynamics 

In this subsection, we will state the results of the study of the traveling waves of lattice 
equations for bistable dynamics in p|] and pQ. In pj], consider a general system of spatially 
discrete reaction diffusion equations for u{t) = {u n (t)} n£ z- 

Un^) = ^ a n,kU n +k(t) + f n (u n (t)), 71 G Z, t > 0, (2.6) 
k 

where the coefficients a n ^ are real numbers and have the following assumptions: 

Al. Periodic medium. There exists a positive integer N such that a n+ N,k = CLn,k and 

/«+*(•) = /»(•) e C 2 (R) for all n, k G Z. 

A2. Existence of ordered, periodic equilibria. There exist ± = {0^} ng z such that 
^an^n+k + /n(0 ± ) = 0, (fi^+N = 4>nAn < 4>ni n e z - After an appropriate change- 

k 

of-variables, the equilibria take the form <ft~ = and <f) + = 1. 

A3. Ellipticity. a n> & > for all k ^ and a n fl = — a n ,fc < 0, n G Z. 

A4. Nondecoupledness. For every integer pair i ^ j, there exist integers ■■■,i m 

rn—l 

such that io = i and i m = j with J^J ai 3i j s+1 _j s > 0. 

s=0 

A5. Finite-range interaction. There exists a positive integer fc such that a„ fc = for 
| A; | > ko and for all n G Z. 



Theorem 2.2. fS'ee /^/j. Assume that and 1 are steady-states and any other N-periodic 
state (p(t) = {(fi n (t)} n( zz with (fi n G (0,1), if it exists, is unstable. Then the problem A2.6\) 
admits a solution (c,w) satisfying vj(-oo) = < to(£) < 1 = w(+oo) for all (el. 

The following theorem is followed by Bates, Chen and Chmaj's results Theorem 1 in 
PP. Consider the following system, 

i ak(u(x + ek) — 2u(x) + u{x — ek)) + f{u) = 0, 



c — (2.7) 



c f u — 

c- 

k>0 

u(±oo) = ±1, 

12 



under the assumptions: (i) / G C (R) has exactly three zeros, —1, q G (—1,1) and 1, with 
/u(±l) > 0; 



5 






(2.8) 



It is well-known that (12.81) has a unique traveling wave solution denoted by (co,0o)- 

Theorem 2.3. (See fl^) Suppose Co 7^ 0. Then there exists a positive constant e* such that 
for every e G (0, e*) ; the problem \2. 7| ) admits a solution (c t ,4> t ) satisfying lim(c e ,0 e ) = 

e— >0 

(c o ,0o) in K x H\R). 

Next, we extend the results in PQ to vector LDEs. 

3 Perturbations of Lattice Differential Equations with 
Bistable Nonlinearity: Persistence of Traveling Waves 

In this section, our goal is to study the persistence of traveling waves of the lattice 
differential equations, 



Ni 

where A is defined as in (12. 5p . that is, Am = —cu'{x) — Aj(x)u(x + rj) with r x = 0, 
and Tj rk,l < j < k < Ni, Ni > 2. The perturbed system of (13. ip is of the form, 



where e > and Bu := Bj(x)u(x+lj) with l x = 0, and lj ^ l k ,l < j < k < N 2 , N 2 > 2. 



We now give the assumptions for the systems of (13.11) and (13.21) . We make the following 
assumption for the nonlinear term: 

(Hl)Fj G C 2 (R), with Fi(0) = 0,Fi(l) = 0,F((0) > and Ff(l) > for i=l,..., N. 

We remark that (HI) is the typical property for the bistable nonlinearity. Assume that 

(H2) There exists a traveling wave solution connecting and 1 for A3. 

We let (c , 0o) be a traveling wave with speed c > for (13. ip . We make the following 
assumption for the perturbed term: 

(113)50 = 0, B\ = and \\B(f) \\ L 2 < 00. 

For simplicity, let A t = A + eB. We may write (I3.2p in 



Am + F(u) = 0, m(+oo) = 1, tt(-oo) = 0; 



(3.1) 



Au + eBu + F(u) = 0, u(+oo) = 1, u(-oo) = 



(3.2) 




A e + F(u) = 0, u(+oo) = 1, u(-oo) = 0. 



(3.3) 
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It is natural to hope that at least for small e, (I3.3P also has a traveling wave solution. 
We prove that the following additional conditions (H4-H5) are sufficient for the existence 
of traveling waves of (13. 3p . We make the following assumption for the operator A e : 

(H4) Fore G (0,1], (A e #, *) < a := ^— wm^Fftl), F{(0)} for fixed a > and 
* G L 2 (R,R N ), and (A e tf, *') = /or * G ^(R,^). 
We remark that according to (HI), a in (H4) is positive. 

Let 7(0o) be a N by N matrix with 7^ = -F/(0o) otherwise 7^ = for i ^ j, and 
Lq := cq4>' — A o + 7(0 O )0 and Lq := — Co0' — Aq0 + 7(0 O )0, where Aq\1/ is the adjoint 
operator of A . We assume that 

(H5) Lq are asymptotically hyperbolic. 

(H5) is equivalent to the assumption that are hyperbolic at ±00. Note that 0o(oo) — 
1 and o (-oo) = 0. Let L+ := c o 0' - A o + 7 (1)0, Lt^ := c o 0' - A o + 7(0)0, 
L~0 := -c o 0'-Aq0 + 7(1)0 and LZ^cj) := — c o 0' - Aq0 + 7(0)0. Then (H5) is equivalent 
to the following: 

(H5) and L^^ are hyperbolic. 

In applications, we may use (H5) instead of (H5) if needed since (H5i) is independent 
of 0o and can more easily get verified. 

Theorem 3.1. Suppose cq 7^ 0. Assume HI — H5. Then there exists a positive constant 
e* such that for every e G (0,e*) ; the problem (13.31) admits a solution (c e , <p e ) satisfying 

lim (c e ,0 e ) = (c o ,0 o ). 

To prove Theorem 13.11 we borrow ideas from [TJ, which are applicable to vector LDEs. 
We remark that the modification is not trivial. We made assumption (H2) for (13.11) instead 
of some specific equation having a traveling wave solution like (12.81) . Theorem 12.21 in Sec- 
tion 2 that Chen, Guo and Wu proved in [4] provides nice candidates for (13. ip satisfying 
(H2). To study the invertibility of some operators, we made assumption (H5) and inves- 
tigate their hyperbolicity, where the Fredholm alternative theory in [12] of Mallet- Paret 
plays an important role. 

Following the ideas as in p] , we let = 0o + ip and formulate the problem as 

L+rf = R(c,4>), (3.4) 

where 

Lt^ = cotf ~ K^P + 7(0o)^ + StP, (3.5) 

R(c, rf>) = (c - c)(0 o + V') + e50 o + 6ij> - iV(0 o , </>), (3.6) 

iV(0 o , V) = F(0o + V) - F(h) - l(<h)il>- (3-7) 
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In some places, we need study the adjoint operator of Lf s , and find its adjoint 

= -co^' - + 7 (0o)^ + Sif>, 

where A* is the adjoint operator of A £ . 

If the inverse of Lf s exists, then let Tip = (L^ h (5 )~ 1 i?(c, ip) and we may rewrite (13.41) as 

Tip = ip. (3.8) 

Then, we can prove the existence of the traveling wave solutions by showing that T is a 
contraction. 

Let ^ o + = 0Vll0olU 2 - 
Lemma 3.1. 

(1) L^ipQ = 0. There exists ip$ G L 2 such that Lqi/jq = with HV'cTlU 2 = 1- Moreover 
ip% G H\R,R N ). 

(2) There exists a positive constant C\, which depends only on F, such that 

for all G H^R, R N ) satisfying (j)±ipQ. 

(3) There exists positive constants C2 such that for 5 G (0,So), 

W\v><<h{\\1>\\v> + 

for all <p G H\R, R N ) where ip = L±(j) + 50 and 5 = ( ^° > . 

Proof. (1) Lqi[)q = follows by differentiating the equation and a direct computation. 
By Theorem 12. 1[ dim(K(L,Q)) = dim(K(LiQ)), and then there exists ipQ G L 2 such that 
LqiPq = with H^cTlU 2 = !• Note that c (ip£)' = A ip+ - ~({(po)ip£ and c (iPq)' = 
-A* ipo + 7(0o)^(7, which imply that 1$ G H 1 (R,R N ). 

(2) If not, there would exist sequences {0„}£° =1 C H\R,R N ) and {ip n }n=i C L 2 (R,R N ) 
such that Lq <p n = ip n , ip n G K^(L Q ), but 

||0n||i/i > 0„|| L 2. 

Without loss of generality, we assume H^nlli? 1 = 1- Thus, we have ||L Q n |U 2 — )• in L 2 
as n —7- 00. As {0 n }^ =1 C H l {R, R N ) is bounded, there exists a subsequence such that 
4>nj ~ ► u m -^ 2 - As are closed, we have L^u = 0, which implies that u G K(L ). On 
the other hand, by the construction, u is in the orthogonal complement of denoted 
by K ± (Lq). Note that K(L^) n K x (Lf) = {0}, then u = 0, which contradicts with 
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(3) The proof follows by arguments similar to those in Lemma 3.1 of pQ. For reader's 
convenience, we provide a proof here. Note that ip = L^<fi + 5<fi and then Lq4> = ip — 5(f). 
Decompose = (4>,ijj^)^ + (j) 1 - with (j) 1 - _L ip^. Thus by (2), 

U L \\m < Ci\\L±^\\v> = Cill^H^ < dfllVlb + 8\\<t>\W). 

Let a = (^,ipo). Note that (i/j - Sfaip^) = (L±<p,^ = (faL^) = 0. Thus 
°\{fa^)\ < ||0 ± IU 2 + ll(^#)|.Thereforea|(0,#)| + ||^||^ < $\(1>,1>*)\+2C 1 (\\il>\\ I ?+ 
5U\\ L 2). Since ||0|| £a < |<0,#}| + ||^|| £a and a < 1, (a - 2(7^) ||0|| L 2 < + 
2Ci||V , IU 2 - Choose 5 < <5o = -^r and C2 = This completes the proof. □ 

For every e, 5 G (0, <5q], we define 




Lemma 3.2. There exists a positive constant Cq such that A ± (5) > 2/Cq for all 5 G (0,<$o]. 

Proof. The proof follows by arguments similar to those in Lemma 3.2 of pQ. By the 
definition of A ± (<5), there exists a sequence {(cj,(f>j)}'jLi in (0, 1) x i/ 1 (M, M^) such that 
lim 6j = 0, = l,for all j, and ijjj := L e . s <fij, satisfies 

j—too " 

lim{||^|| L2 + =A ± («5). 

There exist functions G H^R,^) and ^ G L 2 (R,R^) such that, as j 00, fa (j) 
in L^ oc (R, R^) and weakly in H l (R,M. N ), ipj -»■ ^ weakly in L 2 (M, R^). By the weak 
lower semi-continuity of the L 2 norm, HV'jlU 2 + 1/^1(^5 0o )l < ^■ ± (^)- Thus, for any test 
function <; G C^°(M, W. N ), (ipj,s) = (L eS fa,s) = (fa,L^ s s). Sending j — > 00, we obtain 
(^) = (faLfiq), for all ? G Cq°(R, R^'). That is, <\> G !ff 1 (R,R JV ') is a weak solution to 
L±0 = V- Consequently, by Lemma O (3), ||0|| £2 < C 2 {\\^\\ L , + f|(V>, < C 2 A±(5). 

It remains to find a positive lower bound of ||0|| £ 2- 

Note that (Lfjfa,fa) = (ipj,fa) and (A e fafa) = 0. We have 

Co||0j-|| L 2 = - (7(0o)0j,0j)- 

The Cauchy inequality then gives ll^^o^looll^jll^ > Co||0j||| 2 — IIV^'IU 2 ) which implies 

nn<h)\\Uh\\b> %\mh-ni>j\\h- 

Choose m > such that F-((j) ) > a for 1 < i < N and |x| > m where a is as in (H4). 
Using (ipj,(f)j) = and (A e \l/, < a (see assumption (H4)), (ipj,fa) = (Lf s fa,fa) and 
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we obtain 



{^Ai) > (7(0o)0„0,) > min{ min {F/(0 O )}} - ||^(0o)|| 

\x\>m l<i<N 

= a||0,lli 2 -(a+lin0o)||oo) / 0?- 




J x\<m 



Therefore 



(a + ||F'(0 o )||oo) / 



II 2 



J \x 



x\<m 



Thus, we see that there exist positive constants C3 and C4, which depend on |cq| > and 



Theorem 3.2. There exists a positive constant Co and a positive function 6q such that for 
every 5 G (0, 5q] and e G (0,e ), A ± (e,5) is a homeomorphism for M^) to L 2 (R, R N ) 



for all ip G L 2 with vp J_ tpjf 

Proof. The proof follows by arguments similar to those in Theorem 2 of pQ. Since A ± (<5) > 
2/C , there exists e (<5) > such that A(e, 5) > 1/C for every e G (0,e (5)]. Now we 
consider the operator Lf s with e G (0,eo(5)]. First of all, Lf s is a bounded operator 
from iy^RjR*) to L 2 (R,R N ). Also, by the definition and the lower bound of A(e,5), 
is a homeomorphism from M^) to its image Lf 5 (H 1 (R, M. N )), and the inverse 

(Lf^)- 1 from Lf/H 1 (R,R N )) C L 2 (R,R*) to ff^R*) is bounded by 1/A(e,5) < 
Co- As a consequence of the boundedness, we see that Lf s (H 1 (R,R N )) is closed in 
L 2 (R,R N ). It remains to show that L^ S (H 1 (R, R N )) = L 2 (R,R N ). Indeed, if this were 
not true, there would exist a non- trivial ip G L 2 (R, R N ) orthogonal to Lf s (H 1 (R,R N )), 
i.e., = (Lf^ip) = {</>, L* s i{>) for all G H X (R,R N ). Thus = in H^R"). 

Consequently, Tc Q ip' = (A e - 7(^0) - <S)V is in L 2 (M, R^) since iff G L 2 (M, R^). 

Thus, ip G R ) and L^ s ip = 0. By the definition and the positivity of A(e, 5), we 

then conclude that ip — 0, and therefore obtain a contradiction. Thus Lf s (H 1 (R,R N )) = 
L 2 (R, R N ), and this completes the proof. 




and 




\\(L%r^\\ m <c \m L 2, 



□ 
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Let c(i/j) be the unique constant such that R(c, i[))A-i/; Q so that 
Define 

We claim that T is a contraction. 

Proof of Theorem lff.il The proof follows by arguments similar to those in Theorem 
1 of [1]. First we define X v := {</> G /f 1 (R,R JV ') : ||0|| H i < »7>. We need to show that T 
is a contraction, that is, there exists a A < 1 such that \\Tipi — Tip 2 \\H 1 < Mlrfi ~ ^2 1 If/ 1 
for all V'ij'02 e -^77 ■ By the above Theorem 13.21 ||T^i — T^H// 1 < Coll-R^i — R^Wl 2 - It 
suffices to show that \\R1p1 — R^ 2 \\l 2 < Mlrfi — V^llff 1 , for some A < 1/Cq. 
Recall that 

m= (e£0 o , ^ -) + 5(tfj, y>Q ) - (jV(0 o , ip),i>o) 

Let 5 := |(0o, ipo}{> 0) and then for ip G X,,, ((fi' ,ipo) + {ip'^o) > $ provided 77 <5. 
Note that 

|iV w (0o,^)l <M V \4>®\ 

and 

\N®{<f>oM -iV«(0o,^)l < Mry^f - 1,2, 

where M = max { sup F/'(s)j. Thus we have 

^ N \s\<i+S 

, m _ , , Wo, ^0") + Sty, gT) ~ Wpj ^)^o") I 

and |c(Vi) -c(V 2 )| < - where ^ = t 2 [\\{eB^)\\ L , + (5 + 77) (<f + M V )}. 

Thus, pVi-i#2|U2 < K 2 Ui-Mm, where X 2 = <H<5 + 77 + ||0 o || L2 ][||( eJ B0 o )|| L 2 + 
(5 + 77)(5 + M77)]. Therefore, - T^ 2 \\ m < C \\R^ - Rfa\\ L 2 < C Q K^ X - f 2 \\ m . 

Since ||(e.B0o)||L 2 — > as e — > 0, and we can freely choose 5 G (0, 5q) with <5 as in Lemma 
13.1( 3) and 77 G (0, 5), by choosing < e* <C 1 and appropriate 5 and 77 we can make K 2 
small enough such that A = C K 2 < 1 for any e G (0, e*). This completes the proof. □ 

Theorem 3.3. Suppose Cq 7^ 0. Assume HI — H5. Suppose that for every e G (0, e*), 
the problem admits a solution (c e , e ) constructed as in Theorem \3.1\ Then, there exists 
a subsequence cf) tk such that the problem admits a solution (c e *,0 e *) satisfying 

(c e *,0 e .) = lim (c efc ,0 e J. 
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Proof. Note that |0 e (£)| < 1 and \<p' e (0\ < ||A|| +e*||B|| + max \F(£)\. By Arzela- Ascoli 
theorem, there exists a subsequence u ek that converges uniformly on bounded set. Recall 
that 

Let c e * = lim c e . This completes the proof. □ 

Remark 3.1. Replacing (co,0o) &y (c e *,0e*) and following the arguments of the proofs 
in Theorems \3.1\ and \3.3\ e* can fre extended further unless either (H4) or (H5) is not 
satisfied. 

4 Applications:Existence of Traveling Waves for 1-D 
LDE with Mixed Type 

Let d = d\ + 4d2. 

4.1 the case with d > (1-D) 

Define 

A £ m := — [di(u(x + e) — 2u(x) + u(x — e)) + ^(wfx + 2e) — 2-a(x) + it(x — 2e))l. 

Consider the following equation, 

'cu 7 - A e w + / a (w) = 0, jeZ 

f a (u) = u(u-a)(u-l), (4.1) 
u{— oo) = 0,m(oo) = 1 . 

As e — > 0, we have , 

' cu' - du" + f a (u) = 0, 

/(u) = u(u-a)(u-l), (4.2) 

^u(-oo) = 0,m(oo) = 1 , 

For d > 0, it is well-known that the equation (14.21) has a unique traveling wave O and 

the speed cq. 

Now we consider the system ( 14. ip . By changing variables, we can make f satisfy (i) in 
Theorem 12.31 Equation (14. ip is a particular case with k = 2 of the equation (12.71) . For 
assumption (ii), ^^afc(l — cos(kz)) > for all z £ [0,2tt] is equivalent to the following 

fc>0 

assumption: 

(Al) di + 4d 2 >0. 

Thus, if we assume (Al), both assumptions in Theorem 12.31 are satisfied, and we have 
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Theorem 4.1. Suppose cq ^ 0. Assume that (Al) holds. Then there exists a positive 
constant e* such that for every e G (0, e*) ; the problem (4-1) admits a solution (c e , e ) 
satisfying lim(c e , <f> e ) = (cq, <f>o) mix H 1 ^). 



Remark 4.1. If both d\ and d 2 are positive, then (Al) are satisfied automatically. One 
of them may be negative in the (Al). 

If (Al) is not satisfied, for d < 0, we will transform our model to a new one which is 
in the framework of perturbation method developed in the previous section. In section 
4.2, we will consider the case with d < but d± dominates d 2 in the sense \di\ 3> | g?2 | - 
In section 4.3, we will deal with the case with d < but d 2 dominates di in the sense 
\d 2 \ > 

4.2 the case with d < but \d±\ ^> 1^1 (2-D Transformation) 

As in the work of Brucal - Hallare and Van Vleck |2] , we will use a 2-D transformation. 
First we write the even and odd nodes of the above equation as x = {xj}j^z N an d 
V = {Vj}jGZ N y respectively, and obtain 

Xk = di(y k - 2x k + y k ^i) + d 2 {x k+1 - 2x k + x fc _i) - f a (x k ), j G Z N 
Vk = d\{x k+1 - 2y k + x k ) + d 2 (y k+1 - 2y k + y k _ t ) - f a (y k ) 
To compute the equilibria, define (x±,y±) by 

lim (xj,yj) = \im (xj,yj) = (x + ,y + ) 

The equilibria satisfy E := {(x, y) G M 2 |?/ = x + f a (x) = —f a {y)}- Let 

xj - x_ _ yj - y_ 



(4.3) 



Vj = — ,Wj — . 

x+ -x- y+- y_ 

Then substituting into (14.31) we obtain 

v k = d e {w k - 2v k + w k _i) + d 2 (v k+1 - 2v k + v k _ x ) - f e (v k ), j G Z N 
w k = d Q (v k+1 - 2w k + v k ) + d 2 {w k+1 - 2w k + w k -i) - fo{w k ), 

where d e = d\ v+ ~ v ~ d D = d\ x+ ~ x ~ and 

fe = (X+ ~ xJ) 2 f ae (Vj), fo = (V+ ~ y-) 2 fa {Wj), 

with 

/"(*-) . f"(v-) . 

1, a n = 1. 



(4.4) 



x+ -x- y + - y- 

By choosing proper x,y such that d e , do > 0. If d 2 = 0, this is the case studied in [2]. 
We remark that the case with d 2 = can be easily extended to the case with d 2 > 0. 
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Let * := 0, w) T . Define 
A ^ : = 



±d e (w(x) -2v(x) +w(x- h)), jeZ N 
-jpd (v(x + h) - 2w(x) + v(x)), 



and 

A f # : 



p-[<i e (w(x) — 2v(x) + w(x — h)) + ed 2 (v (x + h) — 2v(x) + v(x — h))], 
p-[<i (f (x + h) — 2w(x) + v(x)) + ed 2 (w(x + h) — 2w(x) + w(x — h))] 

If g(x) = (gi(x), g 2 (x)) T G H 1 ^, M?), we denote g' = (g' 1 ,g' 2 ) T . Now consider 

c^'- A e ^ + F(^) =0, 



and 



(4.5) 

#(-oo) = 0, tf(oo) = 1 , 

c*' - AqVI/ + = 0, 

tf(-oo) = 0,*(oo) = 1 , 
where d e = df^fr, d = d x ^ and = (f e (v), f a (w)) T . 

We can pick those equilibria (x ± ,y ± ) such that after the transformation, any other 
2-periodic state </>(t) = {0 n (^)}ngz with n G (0, 1), if it exists, is unstable. By Theorem 
12.21 there exists a traveling wave solution (cq, o ) for (14. 6ft . To study the system ( 14. 6ft . we 
will apply the perturbation arguments in Section 3.1. We assume 

(Bl)<A e *,*> < a := -^mi n {^(l),^(l),/^(0),^(0)} ; for a > and * G L 2 (M,M 2 ). 

1 + 0" 

Throughout this section, we assume that h G (0, h*) such that (v(x + h) —2v(x)+v(x — 
h),v) < and (w(x + h) — 2w(x) + w(x — h),w) < 0. Note that (A^, is decreasing 
with respect to e for d 2 < 0. It is easy to get the following lemma. 

Lemma 4.1. If (Bl) is satisfied for e — 1, then (Bl) is satisfied for e G (0, 1]. 

Thus, for (Bl), we only need to check for e = 1. 

Let (j) = <po + ip. Following PQ, we formulate the problem as 

= R(c,i/>), 

where 

Lip = c il}' - A e iP + 7 (0o)^ + 
R(c,iP) = (c o -c)(0[ ) + ^ / ) + (A e - A )<po + 5ip - N(<j) ,iP), 
N(<f> , 1>) = F(0 O + ip) - F(<f> ) - T (0o)V. 
In some places, we need study the adjoint operator L* of L, and find its adjoint 

L*ip = -c ip' - A* t i/j + 7 (0 O )V + Sip, 
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where 
A** := 



■jp[d w(x) — 2d e v(x) + d Q w(x — h) + ed 2 (v(x + h) — 2v(x) + v(x — h))], 
■jp[d e v(x + h) — 2d Q w(x) + d e v(x) + ed 2 (w(x + h) — 2w(x) + w(x — h))]. 



Lemma 4.2. For any <p e TT^R 2 ), (A e 0,0') = 0. 

Proof. The proof follows from the direct computation. □ 

To investigate the assumption (H5), we let L e <p := c e (f)' — A e + 7(0 e )0 and 7*0 : = 
— c e (f)' — A*(f) + 7(0 e )0, where c e ^ 0, e (oo) = 1 and 4> e (— oo) = 0. We assume that 
(B2) 7 e (±oo) and 7*(±oo) are hyperbolic. 

We rewrite operators 7 e (±oo)0 := c t (p'— A e 0+7 ± and 7*(±oo)0 := — c e 0'— A*0+7 ± 0, 
where 7 ± = (f 7 ± ) for i=l,2 with 7l + = /:(1),7 2 + = #(1) and 7 f = ^(0),^ = #(0). 

To determine whether 7 e (±oo) or 7*(±oo) are hyperbolic, we need to determine if 
T M±oo) (i0) or T L . (±oo) (i0) ± for 9 G R, where T Le(±oo) (s) := det(s7 - T% =1 A+e ri9 ) 

and T L * (±oo) (a) := detfaJ-E^^V*-) with = ( ° f ) , A\ = ( ^ 7l ) , = 

(% o)> and = (^3 ) T > A 2 = (A 2 ) T ,A^ = (A~l) T , where (A+) T is the transpose of A] 
for i = l,2, 3. First we consider the operator 7 e (±oo). Without loss of generality, we set 
h — 1 and s = i0, we have 

s - [ed 2 ( e -* _ 2 - e s ) - 2d e - 7±] -d e (l + e~ s ) , 

L £ (±oo)lsj - | ^ _ do[l + eS) s _ [ed ^ e _ s _ 2 _ &s) _ ^ _ ^) | 

_ .,s - [2ed 2 (cos9 - 1) - 2d e - 7 ^] -d e (l + e~ s ) v. 

" ^ -d (l + e s ) s- [2ed 2 {cos6 - 1) - 2d - 1 

= -9 2 - i[4ed 2 (cos9 - 1) - 2d e - 7 ± - 2d - 7 ±]0 

+ [2ed 2 (cos0 - 1) - 2d e - 7 f][2ed2 (cos0 - 1) - 2d - 7^] - 2d e d (l + cos9). 

The imaginary part Im(T L ^ ±OQ )(s)) = [Aed 2 (cos9 — 1) — 2d e — 7 f~ — 2d G — 7 2 t ]#. 

The real part Re{T Le{±oo) (s)) = -9 2 + [2ed 2 {cos9 - 1) - 2d e - jf}[2ed 2 (cos9 - 1) - 
2d G — 7^] — 2d e d (l + cos0). We have the following lemma. 

Lemma 4.3. 7/0^0, T Le(±oo) (i0) ^ 0. 

Proof. Let s = i6. If 7m(T Le(±oo) (s)) = and 9 ^ 0, then Aed 2 (cos9-l) -2d e -^f-2d - 
-yi = 0. Thus, 2ed 2 (cos9 - 1) = d e + d + l/2 7 f + 1/2 T |=]. Plugging into i?e(T M±oo) (s)), 
we have i?e(T LE(±oo) (s)) = -0 2 + [d - d e - l/2 7 ± + 1/2 7 ^R - d + l/2 7 f - l/2 7 ±] - 
2d e d G (l + cos9). Since d ci e = d 2 , we have 72e(T i£ ( ±00 )(s)) = — 9 2 — [d — d e — l/2 7 ^ + 
l/2 72 ± ] 2 - 2d 2 (l + cos0) < 0. □ 

Note that Tl € (±o6){s) = ^ l*(±od)(s) since they are symmetric with respect to d Q ,d e . 
Thus, we have the following corollary. 
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Corollary 4.1. The following are equivalent: 

(1) L e (±oo) are hyperbolic; 

(2) L*(±oo) are hyperbolic; 

(3) T Le(±oo) (0) ^ orT i?(±oo) (0)^0. 

We remark that for some typical bistable nonlinearity, we do have T£ e (±oo)(0) ^ or 
T^(±oo)(0) ^ 0. 

Lemma 4.4. If 7^7* > 0, T Le(±oo) (0) ^ and T L * (±oo) (0) ^ 0. 



Proof. We have 



T Le(±O o)(0) = K _ 2do 2do + 7 ±)l 

= 2d e jf + 2^7? + 7i ± 7 2 ± 
> 0. 



□ 



Recall that in !4.6[ f e (u) = (x + —xJ) 2 u{u—a e ){u—\) and f (u) = (y + —y^) 2 u(u—a )(u— 
1). As long e (0, 1), it is easy to verify that f' e {l)f' {l) > and /e(0)#(0) > 0. 

Thus, we have 7f~7^ > and then (B2) or (H5) is satisfied. 

Let c(/0) be the unique constant such that R(c, ip)-LipQ so that 

m= (A e 0o - Aq0 o , jV) + y> -) - (iV(0 o , VQ, y>Q ) 

Define 

T^ = L- 1 i?(c(^),^). 

We claim that T is a contraction. 

Theorem 4.2. Suppose Co 7^ 0. Assume that Bl with e — 1. JTien i/ie problem d4.6[ ) 
admits a traveling wave solution. 

Proof. It suffices to verify the assumptions (H1-H5) in Theorem I3.ll Note that 

fe = (x+ - Z_) 2 /a e 0j), f = (y+ - y^) 2 f ao {wj). 

As long as we choose those f e , f D with bistable nonlinearities, (HI) is satisfied. By the 
choices of equilibria, d e and d D are positive. By Theorem 12.21 (H2) is satisfied. Let 
B = Ai - A . It is obvious that Bl = B0 = 0. Note that 

\\B(j) e \\l 2 = (d 2 f I [v(x + h)-2v(x)+v(x~h)] 2 dx+(d 2 ) 2 / [w(x+h)-2w(x)+w(x~h)] 2 dx 
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Since v(oo) = 1 and v(— oo) = 0, there exists M > such that L X \ >M [ V {% + h) — 2v(x) + 
v(x — h)] 2 dx < 1 and L X \ >M [ W { X + h) — 2w(x) + w(x — h)] 2 dx < 1. Thus, ||50 e ||^ 2 < oo 
and (H3) is satisfied. By Lemma [4.11 as long as (H4) is satisfied for e = 1, (H4) is true 
for all e G (0, 1]. Note that the assumption (H5) is automatically satisfied. By Remark 



13. H e can be extended to 1. This completes the proof. 



□ 



Up to now we studied the cases under d > and d < with \d\\ ^> (g^I- Next we study 
the case d < with \d 2 \ ^> \d\\. 



4.3 (4-D Transformation) the case with d < but \d2 | ^> \d\\ 

If | c?2 1 ^> |c?x|> 2-D transformation may not work because (H4) may not be satisfied. 
Instead of writing the nodes in ordered pairs, we write the nodes of the equation as 
w = {wj}j e %N, x = {xj}j e %N, y = {yj}j£z N an d z — { z j}j&z N i respectively. We obtain 
'w k = di(z k -i - 2w k + x k ) + d 2 {y k -i - 2w k + y k ) - f a {w k ), 
x k = di(w k - 2x k + y k ) + d 2 {z k - 1 - 2x k + z k ) - f a (x k ), 
y k = di(x k - 2y k + z k ) + d 2 (w k - 2y k + w k+1 ) - f a (y k ), 
yz k = dx(y k - 2z k + w k+1 ) + d 2 (x k - 2z k + x k+1 ) - f a {z k ), k G Z N . 
To compute the equilibria, define {w±,x±,y±, z±) by 



(4.7) 



lim («>■'>//;• %) = (w_,a;_,y_,z_), lim (wj, Xj, y jt Zj) 

J— >— oo j— >oo 

The equilibria (w, x, y, z) G M 4 satisfies 

' d x {z -2w + x) + 2d 2 (y - w) = f a {w), 
di(w - 2x + y) + 2d 2 (z - x) = f a (x), 
d x {x -2y + z)+ 2d 2 (w - y) = f a (y), 
K d x (y -2z + w) + 2d 2 (x - z) = f a (z). 



[w+,x + ,y + ,z + ; 



(4i 



Let 



W ; 



Wj — W„ 



, Xj 



yj - y- 



"J Z 3 



Zj - Z_ 



w+ — W- x + — X- U+ — V- z+ — Z- 

Then substituting into (14.71) . for simplicity, we discard the hats for w,x,y,z and obtain 
' w k = a u z k ^t - b u w k + b 12 x k + a 13 y k ^ + b 13 y k - fi(w k ), 
ik = b 21 w k - b 22 x k + b 23 y k + a 2i z k -i + b 24 z k - f 2 (x k ), 
Vk = b 32 x k - b 33 y k + b u z k + b 31 w k + c 31 w k+1 - f 3 (y k ), 
{z k = b 43 y k - b 44 z k + c 41 w k+1 + b 42 x k + c 42 x k+ i - f 4 (z k ), k G Z N , 



(4.9) 



where Ai 



A 2 = (bij) and A 3 = (c^j) are 4 by 4 matrices given by 



/0 d 2 ^^ 



A, 



W+—W. 





\0 









d-i 



d 2 ^- 

^ X+—X- 





/ 
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A, 



J w 



hi 

W+—W- 



X+—X- 



d *±^ 



y+-y- 



X+—X- 
W-\-—W- 

-b 22 

d X ^^ 

y+-y- 



' W+—W- 

d x ^± 

1 X+—X- 

~b 3 3 



d 2 ^ 

Z X + 



A, 






d 2 x+ ~ 


~ x - dt y+ - 

-z- 1 


-v- 








-Z- 






( 


























d w+-w- 

y+-y- 















d 2 x+ - x - 





0/ 



y+-y- 
-6 44 



with bu (i=l,2,3,4) is such that A 1 + A 2 + A 3 = 0. Note that {A x + A 2 + A 3 )l = 0, 
since 1 is an equilibrium, we have = 0. Obviously, /i(0) = 0. Thus, is of form 
fi(0 = ki£,{Z, - - 1) for some k u a G R. 

For simplicity, let x(£) = (w(£), x(£),y{£), z(£,)) T and F{x 
then we consider the following system of equations in 

3 

cx'(0 - MZMZ + rj) + F(x) = 0, x(-oo; 
i=i 

where 7*1 = —h, r 2 = and r\ = h. 



(h(w)J 2 (x),f 3 (y),U(z)) T , 
0,x(oo) = l, (4.10) 



Choose the equilibria with w + — w~ > 0, x + — X- > 0, y + — y- < 0, and z + — z_ < 
(or with w + — u>_ < 0,x + — x_ < 0,y + — y- > 0, and z + — z^ > 0). Without loss of 
generality, assume d\ < 0. We rewrite the system of equations (14.101) as the following: 

3 3 

cx\0-Y, A ^MZ + rj) + J2 B ^MZ + rj) + F(x) = 0,x(-oo) = 0,x(oo) = 1, (4.11) 
j=i 3=1 
where r\ = —h, r 2 = and r% = h, 



A, 



/0 d 2 ^^ 





\0 



A- 2 



' -^1 


d 2 ^^ 

y+-y- 




V 



.4, 





-b 22 


d, x -^ 



Z+—Z- 



( 



\d 






d ™+^ 

• y+-y- 

W+—W- 
Z+—Z — 



di 



d 2 ^ 



Z+—Z- 
W+—W- 



\ 



-x + - 





J 



' W+—W- 







z +- 



-b 



di 





d 2 x+ - 



33 

y+-y- 



d< 

* X+—X 



Z+—Z- 



y±-y- 
-644 



Z+— Z- 



0\ 


0/ 
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and Bi = B 3 



(- 



B 2 = di 



W-^—W- 
W-\-—W — 



W+—W — 

y+-y- w+—w- y+-y- 





V 



X + — X— X+—X- 
X+—X — 



y+-y- 




X+— X — 
X+—X- 

y+-y- 




0\ 






with bu (i=l,2,3,4) is such that A 1 + A 2 + A 3 + B 2 = 0. 
Consider the reference system, 
3 

cx'(0 - MtXt + r j) + F ( x ) = °> *(-°°) = 0. *(°°) = I 

where r\ = —h, r 2 = and r\ = h. 



(4.12) 



We can pick those equilibria (w , x , y , z ) such that after the transformation to 
and 1, any other 4-periodic state (pit) = {0 n (£)}nez with (fi n G (0,1), if it exists, is 
unstable. In this paper, we only consider the cases having bistable dynamics after the 
transformation. By Theorem 12.21 there exists a traveling wave solution (co, 4>o) for (14.121) . 

Denote Ax := Yf j= i Aj{g)x(£ + rj). Let Bx := Y? j= i B j(O x (C + r i) and A e = A + eB. 
We consider the following: 

cx' - A e x + F(x) = 0, x(-oo) = 0, x(oo) = 1. (4.13) 

We assume 

(Dl) For e G (0, l],(A e tf, < a := min {f-(l), //(0)} /or a > and * G 

1 + cr i<«<4 

L 2 (M, M 4 ). 

Let Lq := c o 0' — A <f) + 7(0o)</> and Lq <p := — c o 0' — Aq0 + j(<fi )(f), where Aq\I/ is the 
adjoint operator of A . We assume that 

(D2) Lq are asymptotically hyperbolic. 

Theorem 4.3. Suppose Cq 7^ 0. Assume (Dl) — (D2). Then there exists a positive 
constant e* such that for every e G (0,e*], the problem (I4.13P admits a solution (c e ,<p e ) 
satisfying 

lim (c e ,<j) c ) = (c o ,0 o )- 

By Remark 13. 11 we can extend e* if (H5) is satisfied. If we can extend e* to 1, then we 
successfully get the existence of traveling wave solution for Equation ( 14. lip . We remark 
that, unlike the previous 2D case, (H5) is not automatically satisfied for this 4D system. 
At the end of each extension, we have to check the assumption (H5). Let L e (±oo)0 := 
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c e 0'- A e + 7 ± and L*(±oo)0 

fori=l,2,3,4with 7 + = //(l),7r 
(D2) L e (±oo) and L* (±00) are hyperbolic. 
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